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Abstract 

The boundary conditions for the ideal MHD equations on a plane dis¬ 
continuity surface are investigated. It is shown that, for a given mass flux 
through a discontinuity, its type depends only on the relation between inclina¬ 
tion angles of a magnetic held. Moreover, the conservation laws on a surface 
of discontinuity allow changing a discontinuity type with gradual (continu¬ 
ous) changes in the conditions of plasma how. Then there are the so-called 
transition solutions that satisfy simultaneously two types of discontinuities. 
We obtain all transition solutions on the basis of the complete system of 
boundary conditions for the MHD equations. We also found the expression 
describing a jump of internal energy of the plasma howing through the dis¬ 
continuity. Firstly, this allows constructing a generalized scheme of possible 
continuous transitions between MHD discontinuities. Secondly, it enables 
the examination of the dependence of plasma heating by plasma density and 
conhguration of the magnetic held near the discontinuity surface, i.e., by 
the type of the MHD discontinuity. It is shown that the best conditions for 
heating are carried out in the vicinity of a reconnecting current layer near 
the areas of reverse currents. The result can be helpful in explaining the 
temperature distributions inside the active regions in the solar corona during 
hares observed by modern space observatories in soft and hard X-rays. 
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1. Introduction 


Space observations of solar flares in soft and hard X-rays with space¬ 
craft Yohkoh and RHESSI at hrst revealed the characteristic strnctnre of 


active reg i ons in the solar atmosphere related to flares (iTsnneta et al.l. Il992 


Lin et al.l. [20021). Primary energy release in a flare occnrs at the top of the 


magnetic held loop strnctnre, the base of which extends below the photo¬ 
sphere. Charged particles, accelerated in a hare, fall along the magnetic 
held lines to the snrface of the Snn and collide with the dense chromo¬ 
spheric plasma. Deceleration of the particles is accompa nied by the hard 


X-ray emission at the loop footpoints (ITsnneta et ah . 19921) . Chromospheric 


plasma heated by the collision rises along the magn etic held lines and pro- 
dnces soft X-ray emission of the loop flTsnnetai. Il9961) . Hard X-ray sonrces at 
(or above) the top of the loop are associated with the thermal plasma heated 


direc t ly in (or from) the r e gion of primary ener gy release flMasuda et al. 
1994 : Petrosian et ah . 2002 : Sni fc Holman . 2003 ). 


The energy sonrce of a hare is a non-potential part of magnetic held in 
the solar corona. The so-called free magnetic energy of interacting mag- 
netic hnxes i s accnmulated in the magnetic helds of coronal electric cnrrents 


(ISvrovatskiil. Il962l: iBrnshlinskii et al.l. 119801) . Disrnption or qnick dissipation 


of snch cnrrent strnctnre can lead to the energy release as kinetic energy of 
accelerated particles or thermal energy of heated plasma. The area of energy 
release is described by the reconnecting cnrrent l ayer models, m ore exactly. 


the model of snper-hot tnrbnlent-cnrrent layers f Somovl. 2013b . Chap. 8) 


The high-speed plasma hows form a system of discontinnons hows near a cnr¬ 
rent layer. It is observed, fo r example, in nnmerica l simnlations of th e mag¬ 


netic reconnection process (Shimizn fc Ugail. l2003l: lUgai et al.l. l2005l: lUgai . 


20081: IZenitani fc T. Mivoshil . 120111) . The MHD discont innit ies are also able 
to convert the energy of the directed motion of the plasma with frozen mag¬ 
netic held to the thermal energy, thns making a farther contribntion to the 
heating of a snper-hot plasma in a hare. 

An abrnpt change (jnmp) of p hysical characteristics of a plasma how oc¬ 
cnrs on the discontinnity snrface f Svrovatskii . 1957 : Anderson, 19631) . It may 
be jnmps of density and velocity of the plasma or jnmps of intensity and direc¬ 
tion of the magnetic held lines. The relations between these jnmps determine 
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the type of discontinuity. Ordinary hydrodynamics allows only two types of 
discontinuous flows: a tangential discontinuity and a shock wave. However, 
there is much variety of possible discontinuity types in magnetohydrodynam¬ 
ics (MHD) owing to the presence of the magnetic held. Furthermore, it is 
possible to transit from one how regime to anoth er with a gradual (continu¬ 
ous') changes in the characteristics of the plasma fjSvrovatskiil . Il956l: iPolovinl . 
196ll) . Then one system of boundary conditions for the MHD equations on 
the discontinuity surface would satisfy once two types of discontinuous hows 
at the moment of transition. This is so-called transition solution. 

We aim to study the possibility of the transitions between diherent types 
of MHD discontinuities and plasma heating by the discontinuous hows. Ini¬ 
tially, in Sect. [21 we describe the s t andard classihcation of MHD discontinu¬ 
ities (e .g., Svrovatskii . 1957 : Priestl. 19821: Goedbloed &: Poedts . 2004 : Somov . 
2013al) with a view to associate it with the amount of the mass how in Sect. [3l 
On this basis, we found all transition solutions for the full system of bound¬ 
ary conditions (Sect. H]) and construct a demonstrative scheme of allowable 
transitions in MHD (Sect. [5]). Then we investigate the possibility of plasma 
heating by diherent types of MHD discontinuities (Sect. [H]). Finally, we dis¬ 
cuss our hndings as possibly applied to calculations of the analytical model 
of the magnetic reconnection in the context of the basic physics of solar hares 
(Sect. [7]). 


2. Boundary conditions 

We will seek a solution of the formulated problem for an MHD discontinu¬ 
ity, i.e., a plasma region where the density, pressure, velocity, and magnetic 
held strength of the medium change abruptly at a distance comparable to the 
particle mean free path. The physical processes inside such a discontinuity 
are deter mined by kinetic phenomena in the plasma , both laminar and turbu- 


? K 

lent ones fiTjOngmirel. Il963t iTideman &: Kralll. Il97ll) . In the approximation of 


dissipative MHD, the internal structure of a discontinuous how is dehned by 
dissipative transport coeffi cients (the viscosity and electric conductivity') and 
the t hermal conductivity (jSirotina fc Svrovatskiil . Il960l: IZeldovich &: Raizerl . 
19671) . However, in the approximation of ideal MHD, the jump has zero 
thickness, i.e., it occurs at some discontinuity surface. 

We will consider a plane discontinuity surface, which is appropriate for 
areas of a sufficiently small size compared to the radius of curvature of the 
discontinuity surface. Let us introduce a Cartesian coordinate system in 
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Figure 1: Changes in the magnetic field B, velocity field v, and plasma density p at the 
shock front x = 0. Since ^ 0 in the coordinate system associated with the current 
layer, the velocity vectors are not parallel to the magnetic field vectors. 


which the observer moves with the discontinuity surface located in the (t/, z), 
plane in the direction of the x axis (Fig. [1]). In the approximation of ideal 
MHD, we neglect the plasma viscosity, thermal conductivity, and electric 
resistivity. The boundary conditions for the MHD equations at the discon¬ 
tinuit y then take the form of the following conservation laws flSvrovatskiil . 


19571 1: 


= ( 1 ) 

{pv^} = Q, ( 2 ) 

{ v^By - VyB^ } = 0 , (3) 

{ - v^B^ } = 0 , (4) 

P '^x^y B^By !• 0, (5) 


P BxB^ ^ 0 , 

B ' 

p + P'^x+^ r“ = 


( 6 ) 

(7) 
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I pVx p) ^ - (V • B) Sa;) I = 0. (8) 

Here, the curly brackets denote the difference between the values of the 
quantity contained within the brackets on both sides of the discontinuity 
plane. For example, Eq. ([1]) implies the continuity of the normal magnetic 
field component: 

{Bx} = Bx 2 — Bj-i = 0 . 

The quantities marked by the subscripts “1” and “2” refer to the side corre¬ 
sponding to the plasma inflow and outflow, respectively. 

In contrast to the boundary conditions in ordinary hydrodynamics, the 
system of boundary conditions ([I])-® does not break up into a set of mu¬ 
tually exclusive groups of equations and, hence, in principle it admits con¬ 
tinuous transitions between different types of discontinuous solutions as the 
plasma flow conditions change continuously. Since a smooth transition be¬ 
tween discontinuities of various types is possible, the local external attributes 
of the flow near the discontinuity plane are taken as a basis for their classifi¬ 
cation: the presence or absence of a mass flux and a magnetic flux through 
the discontinuity, the density continuity or jump. 

In the presence of transition solutions, the classification of discontinuities 
in MHD can only be relative. Indeed, a discontinuity of a given type can 
continuously pass into a discontinuity of another type as the plasma inflow 
and magnetic field parameters change gradually. As it will be shown in the 
Sect. [7] the type of discontinuity can change when passing to a different point 
of the discontinuity surface. In any case, since a smooth transition is possible 
between discontinuities of various types, the local external signatures of the 
flow near the discontinuity plane are taken as a basis for their classification: 
the presence or absence of velocity, Vx, and the magnetic field, Bx, compo¬ 
nents perpendicular to the plane (i.e., normal), the continuity or jump in 
density p. With respect to these signatures, the energy conservation law (jS]) 
is an additional condition: at the magnetic field strength, the velocity field, 
and the density jump found, Eq. ([8]) defines the jump in pressure p. Thus, 
considering our objective of identifying the discontinuities near the recon¬ 
necting current layer in mind, we can restrict our analysis to the remaining 
seven equations: dlD-O- 

It is necessary to establish conditions under which certain types of dis¬ 
continuous are formed. Let’s rotate the coordinate system about the x axis 
in order to make Vzi = 0. The substitution of Eqs. ([I])-([2]) in ([6]) will then 
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yield the equation 


{Bz} = pv^ Vz2 ■ 
47r 


(9) 


Here are some simple particular solutions of this equation: 

1. If Bx = 0, Ua; = 0, i.e., there are no magnetic flux and plasma flow 
through the discontinuity, then {vy}, {By}, {Bz} and {^ 22 } are arbi¬ 
trary quantities as we see from Eqs. ([3]), (E]) and (jH])- The Eq. ([7]) 
implies that the pressure and the magnetic held strength are related 
by the continuity of the total pressure 



This solution corresponds to the classical tangential discontinuity. 

2. If = 0, Ua; 7 ^ 0 (then from Eq. [9]nz = 0), the magnetic held is parallel 
to the discontinuity surface and increases with the compression of the 
plasma due to a freezing condition, as we see from Eqs. ([2]), (jl]). 



This is a perpendicular shock, that is well-known in MHD. 

3. If Ha; 7 ^ 0, Ua; = 0 (then from Eq. [9] {B^} = 0), Vy, and By are 
continuous from Eqs. ((3]), (E]). The direction of the magnetic held does 
not change, however, the value of the density may change. This is a 
contact discontinuity. 

4. If Bx 7 ^ 0, Ua; 7 ^ 0, then substitution = 0 in Eq. (E]) gives {B^} = 0. 
As a result, Eq. (|1]) is transformed into 


Bz{vx} = 0 . 


It allows two diherent solutions: 

(a) At the hrst, we consider the solution {v^} = 0. The substitution 
of this condition in Eq. (E]) gives Alfven wave ({p} = 0). Then 
Eq. ([7]) rewriting as {By} = 0. 

(b) Solution Ha, = 0 leads ns to the two-dimensional picture of the 
discontinuity: the velocity and the magnetic held lie in the same 
plane orthogonal to the discontinuity plane. 
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Thus, the presence or absence of the plasma overflow and penetration of 
the magnetic held through the discontinuity surface allows identifying all 
types of discontin u ities listed abo ve; for details of their properties see, e.g., 
Svrovats^ f 19561) : Somov f 2013al) . The description of many discontinuities 


is simplihed in the coordinate system, often called the de Hoffmann-Teller 
system. Such coordinate system moving along the discontinuity plane with 
velocity 


u. 


= V 


— B, 


Then the vector B becomes parallel to the vector v. In particular, Alfven 
discontinuity describes as the rotation of the magnetic held vector around the 
X axis and all oblique discontinuous (with the magnetic held inclined to the 
discontinuity plane) become two-dimensional in this case. The de Hohmann- 
Teller system is inapplicable to tangential discontinuity and perpendicular 
shock wave due to the limitation 7 ^ 0 , and we will not use it. 

However, taking into account that all discontinuities apart a tangential 
and an Alfven discontinuity and a perpendicular shock can be reduced to 
a two-dimensional form, let us begin consideration of the boundary condi¬ 
tions ([I])-(IH]) with the two-dimensional case 4(b). We will study only the 
classihcation features discontinuous hows to Sect. E Therefore, we will not 
use Eq. (jH]). Two-dimensional discontinuities will be described by the hve 
boundary conditions: 


{H„} = 0, {pv^] = {), i^pv^Vy -Q, (10) 

{v^By - VyB^} = 0 , Ipv^ + p + -^ I = 0 . 

In the next section we will study the change of the magnetic held on the 
discontinuity surface based on the system (fTOj) . 


3. Inclination angles 

We will write the system of Eqs. (ITU]) in the linear form with respect to 
the variables {'Vy}-, and {By}\ 

B 

{ua;} - m{r} = 0 , m{vy} - ^ {By} = 0 , 

47r 
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( 11 ) 


By{vx] - Bx;{vy) + mr{By] = 0 , 

m{vx] + {p] + ^{By} = 0 , 

where we introduced new variables r = 1/p and m = pv^- The mean of the 
two quantities is denoted by the tilde r = (ri + r 2 )/ 2 . 

For existing nontrivial solutions of the linear system of Eqs. ca, the 
determinant composed of its coefficients must be equal to zero: 


-1 

0 

m 

0 

0 

m 

0 

-Bx/^t^ 

By 

-Bx 

0 

mf 

m 

0 

{p}/{r} 

By /An 


Let us expand the determinant: 


{P} (B, 


{r} ydvr 




+ m' 


Bl 


+ 




dvr dvr 


m'^r 


0 . 


The latter equation imposes constraints on the admissible values of the mass 
flux m: 


{p} m^ — B^/4:7Tr 
{r} ^2 _ _j_ ^ 2 '^ / 4:71 r 


( 12 ) 


We can not describe the internal structure of the discontinuity in terms 
of ideal MHD. Therefore, the discontinuity appears to the surface of zero 
thickness. However, the irreversible processes in the shock waves leads to 
the entropy increase. It is essential that the amount of the entropy change 
does not depend on the specific dissipation mechanism and its jump is de- 
termined by the mass, m omentum and energy conservation laws (see, e.g., 
Zeldovich fc Raizer . 19671) . The entropy increase according to the Zemplen 
theorem forbids the rarefaction shock waves. Thus, 


w = 



Pi - p2 
P 1 P 2 


< 0 . 


The quantity can not be negative. Therefore, two inequalities must hold: 
either 

dvrr ’ 


< 


(13) 
















or 


m > 


Anr 

The solution of the system of Eqs. (HD is 


(14) 


{v^} = Cm ( ^ - mV ) , 


(15) 


K} = C'm^^, 

W = v(§-mV) , 

{By} =Cm^By, 

We hnd the constant C by substituting the expressions for {vx} and {By} 
for the last equation of the system flTT]) : 


^ {p] B^^+B^y ,, 

C =-T ^- - — m r 


-1 


m‘ 


dvr 


Consider the last two equations of the system ([T5]). Eliminating the con¬ 
stant C from them, we will obtain the relationship between the tangential 
magnetic held components 




m^{r} ~ 

(5J/dvr — m^r) ^ 


Expanding the relations 


{By} — By 2 - Byl 


and 

By = 2^By2 + Byl) , 

we then have 

2 {B^/ An — m'^r) + m'^{r} 
2 {B^/ An — m'^f) — m‘^{r} 


( 16 ) 
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Let us divide both parts of Eq. fllbp by Bj. to obtain the relation of the 
angles between the magnetic held vector and the normal to the discontinuity 
surface on both its sides: 


tan 62 


2 ( / Att — r) + m^{r} 
2 ( B^ / Att — m‘^f) —m2{r} 


Here, ta.n 6 = By/B^, r = 1/p, and the tilde marks the mean values of the 
quantities, r = (ri + r 2 )/ 2. Rewriting this formula by expanding the jumps 
{r} and the means r, we obtain 


tan 02 


• Attvi/B^ — 1 

• A7ir2/B^ — 1 


tan 01 . 


Denote = B^/Atiti and = B^/A 7 ir 2 ; as will be shown below, ruofr 
and rUon correspond to the mass hux through the switch-off and switch-on 
shocks. Note that rrios < ’^on, because according to Zemplens theorem, 
^2 < The formula for the held inclination angles takes form 

(17) 


Let us denote = B^/Anr and m]_ = By I dvrr. Since ^2 < r < ri, 
ruoff < < rrion- Then we will write conditions flT^ and flTT|) as: 


9 2 

m < m, , 

A ’ 


(18) 


2 ^ 2 I 2 

m > + m_L 


(19) 


Note that ruofr, uia and rrion are not independent but are related by the 
equation 


m 


2 

A 




^off 


m, 


2 

oflf 


( 20 ) 


This can be easily verihed by expanding the mean r in the dehnition of 
Based on these results we will consider the properties of discontinuous 
hows; more specihcally, we will establish the possible transitions between 
them. 
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4. Transition solutions 


We begin to seek transition solutions with a search for the conditions of 
possible transitions between various types of two-dimensional MHD flows 
{vz = 0, -02 = 0), i.e., flows for which the velocity field and the magnetic 
field lie in the (x, y) plane. We call such discontinuous flows plane or two- 
dimensional ones. Then, we will hnd the transition solutions that correspond 
to them and establish the form of the solutions that are the transition ones 
to three-dimensional discontinuous flows. 

Equation flT7|l along with conditions flT5]l - flT^ describes the dependence 
of the magnetic held inclination angles on the mass hux through the discon¬ 
tinuity. This dependence can be specihed either by the two parameters ruofr 
and ruon, or, for example, by the quantities pi and {p}. Since we are in¬ 
terested in the classihcation attributes of discontinuities (i.e., the qualitative 
changes of the relation between the angles 9i and 6*2 when varying m^), for 
the time being, we will consider Eq. fll7p without any specihc application 
to certain physical conditions in the plasma. We will choose the parameters 
from clarity considerations. Let and be related as 3 : 4 : 6 . 

We will measure the square of the mass hux in units of m^/4; then = 3 
and = 6 . The dependence 6 ^ 2 (m^) is shown in Fig. |2] or three values 
of the angle 61. The corresponding curves behave identically. Firstly, they 
intersect at one point at = m^g, 62 = 0 here. Secondly, 62 ^ — 9 i when 
—?• = 4. for each curve. Thirdly, they all have a region that does not 

satisfy conditions (flSj) and ([19]), located near 

Let us separate out the regions in Fig. |2]each of which is characterized by 
its own behavior of the dependence of 6*2 on m^. In region I (0 < < m^g), 

the tangential component By2 of the magnetic held vector B 2 behind the 
discontinuity surface decreases with increasing m^. In this case, 0 < 02 < 
01 , i.e., when passing through the discontinuity surface, the tangential held 
component weakens but remains positive. At m? = m^g when crossing the 
discontinuity plane, By2 becomes zero. In region II (m^g < rn? < m^) By2 
is negative and increases in magnitude, but —0i < 02 < 0. In region III 
(m^ < w? < m^^), just as in region II, By changes its sign when crossing the 
discontinuity plane. Now, however. By increases in magnitude (02 < —0i), 
remaining negative. Finally, in region IV {m? > the magnetic held is 
amplihed (02 > 0 i) with its tangential component retaining the positive sign. 

Thus, we have shown precisely how the behavior of the relation between 
the magnetic held inclination angles and, consequently, the type of MHD 
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Figure 2: Inclination angle of the magnetic field 62 behind the discontinuity plane versus 
the square of the mass flux rn? at various angles 0i . The cases of = 1°, = 25°, 

9i = 45° are designated by the solid, dashed, and dotted lines, respectively. 


discontinuity changes with the increasing mass flux. The regions I and II 
correspond to the slow shocks that, respectively, do not reverse (Sl) and 
reverse (Si) the tangential held component. The regions III and IV corre¬ 
spond to the trans-Alfven (Tr) and fast (5+) shocks, respectively. Shock 
waves c orresponding to th e regions II and III is often called intermediate 
shocks f Sterck et ah . 1998h . The plasma velocity is super-Alfven upstream, 
and sub-Alfven downstream at the intermediate shocks while the slow shocks 
is everywhere sub-Alfven and fast shocks is everywhere super-Alfven. We use 
the mass flux as a classifier rather than the plasma velocity. This classifier 
divides the intermediate shocks into two types which have the different prop¬ 
erties. Therefore, we will use different names for the discontinuities in the 
regions II and III. The transition solutions for the discontinuities correspond¬ 
ing to the adjacent regions are realized at the mass flux demarcating these 
regions. 

Consider the behavior of the function 62 {m‘^, 6 i) near the boundary of 
regions II and III, where = m^. The domain of definition of the function 
6 * 2 (m^, 6 * 1 ) to the left and the right of is specified by conditions (1T8|1 . and 
OSli. respectively. In view of Eqs. lira cLiici (I20p . tcLii 02 —^ —tcLii 01 wliGn 


12 
































i.e., By —)■ 0. Inequality 0191) transforms to in this 

case. Therefore, the function 02(m^, 6 i) in region II and region III is dehned 
near m^. However, the right part of the condition also increases with ([19]). 
The equality + m]_ is established at some value of in region III 

and the strongest trans-Alfven shock (increasing the magnetic energy to the 
greatest extent) takes place. As the mass flux increases further, can not 

satisfy conditions fflS]) and flT^ until again becomes equal to + m]_. 
This occurs in region IV, where the strongest fast shock is observed. 

Let us derive the equation of the curve bounding the function 02 0i), 

and, hence, the strongest (for given plasma parameters), fast and trans- 
Alfven shocks. Setting equal to the right part of condition fll4p . we hnd 

Byl = ±2 4:711^171“^ - Bl - By2 , 

where the plus and minus correspond to regions IV and III, respectively. 
Dividing the derived equation by B^, we have 


tan 01 


±2 




tan 02 


( 21 ) 


Substituting Eq. fl^ in fllTl) . we obtain the equation of the sought-for curve 


tan 02 


m? I — 1 
j m\ — 1 


( 22 ) 


The corresponding curves are represented in Fig. [2] by the thin lines. 

Now, we will seek transition solutions in order of the increasing mass flux 
m, starting from m = 0. Let’s consider the transition between the contact 
discontinuity {C) at = 0 and the slow shock in region 1 . The boundary 
conditions for two-dimensional discontinuities follow from Eqs. ([I|)-(I7]) when 
Vz = 0 and Bz = 0 are substituted: 

r H ^ I 

{B^} = 0, {pu^} = 0, }p + pv^ + -^>=0, (23) 


{ '^xBy VyBx 0 , ^ ^ B^By j 0 . 

The solution of these equations in region I presented in Fig. [2] corresponds to 
the slow shock {Sl) that does not change the sign of the tangential magnetic 
held component. This can be easily verihed at small m^. It follows from 
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Eq. fll7l) which, of course, remains applicable under conditions of the two- 
dimensional discontinuous flows that 

/ \ 

tan 6*2 ~ 1- —I- TT tan 6 *i , 

V "^on/ 

i.e., 0 < 6*2 < 6 * 1 , as it must be in the slow shock. Moreover, when m? —)■ 0, 
i.e., Vx —t 0 , it follows from the system ( 12 ^ that {vy} —)■ 0 and { ^ 0 

which is the only limiting case for the slow shock. 

It remains to show that when Vx —)■ 0, conditions (12^ transform to the 
boundary conditions at the contact discontinuity. Indeed, the substitution 
of Ua; = 0 in fl23|) gives 

{ft} = {t'»} = {B»} = {p}=0, (24) 

At the contact discontinuity, the jump in density {p} is nonzero. Other¬ 
wise, all quantities remain continuous. Thus, solution fl24p simultaneously 
describes the slow shock in the limit Vx 0 ,and the contact discontinuity, 
i.e., it is the corresponding transition solution. 

When crossing the boundary of regions I and II, the tangential magnetic 
field component changes its sign. The slow shock (Sl) that does not reverse 
the tangential field component turns into the reversing slow shock (Si). 
The transition solution is realized on the boundary of the regions, when 
62 = 0. The substitution of By2 = 0 in fl23|l gives the corresponding transition 
solution: 

B ^ 

{B,,} = 0, {pu^} = 0, [p + pv^} = -^, (25) 

{ '^y } ^XlByl , PVx ■{ Vy ]■ BxByl . 

Eliminating {vy} from the last two equations, we find 

=Pi^xV4vr = m 2 g, (26) 

which was to be proved. This mass flux corresponds to the switch-off 
shock (S'off): the tangential field component disappears behind the disconti¬ 
nuity plane. This occurs irrespective of the angle 6*1 and corresponds to the 
intersection of the curves at m? = in Fig. [2J 


14 





The reversal of the tangential held component at the boundary of regions 
II and III can be a special case of the three-dimensional Alfven discontinu¬ 
ity (A). Since there is no density jump at the Alfven discontinuity, let us 
substitute {p} = 0 in fl23l) : 

{p} = {B,} = {t>,}=0, |p+A|=o, (27) 

Bx { '^y }■ '^X ■{ By } , P Vx ■{ Vy ]■ — Bx ■{ By j" . 

4:71 

If {By} = 0, then all quantities are continuous and there is no discontinuity. 
Let {By} 7 ^ 0. Eliminating the ratio {By}/{vy} from the last two equations, 
we obtain 

= pB"^/ dvr = . (28) 

When substituting {p} = 0 in ([I])-([7]), we hnd the boundary conditions at 
the Alfven discontinuity: 

{p} = {B.} = {t>«}=0. |p+A^tAAo, (29) 

Bx { '^y } '^X { By } , P Vx ■{ Vy } — Bx ■{ By j- , 

47r 

Bx ■{ 'Vz }■ '^x ■{ Bz: } , p Vx ■{ V z ]■ ~ Bx "{ Bz ]■ • 

47r 

Comparison of the systems fl27D and shows that the boundary conditions 
fl27|l describe the transition discontinuity between the slow shock in the limit 
{p} —)■ 0 and the Alfven how at Vz = 0 and Bz = 0. The discontinuity is 
a special case of the Alfven discontinuity that reverses the tangential mag¬ 
netic held component. Trans-Alfven discontinuities reverse and enhance the 
tangential held component. They occupy region III and are adjacent to the 
Alfven mass hux (|28|) on the right side. The conditions for the transition to 
the Alfven discontinuity are identical to (HT} . 

There can be no how near the boundary of regions III and IV in some 
range of mass huxes. For this reason, the transition between the trans-Alfven 
and fast shock is forbidden. The range narrows as the initial inclination 
angle of the magnetic held decreases to = 0 (Fig. [3]). The strongest fast 
shock takes place at the minimum possible mass hux that is admissible by 
condition ([19]). As the mass hux increases, the tangent of the held inclination 
angle behind the discontinuity plane decreases, asymptotically approaching 
tan 02 = {pi/Pi) tan 01 . 
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Figure 3: Schematic behavior of the dependence at = 0 (a), — 0 (b), and 

M = o(c). 

5. The scheme of transitions 

Varying pi, {p} and leads to contraction or extension of the curves 
presented in Fig. [2] along the coordinate axes without any change of their 
overall structure. For zero 6 i, and {p} the behavior of dependence 02 
is shown in Fig. [3l In view of Eq. flT7|) when 0i —)■ 0, the angle 62 also 
approaches zero at almost all values of (the case of will be 

considered separately). If Byi = 0, then By 2 = 0 (Fig. [3]a). In this case, the 
boundary conditions for the two-dimensional discontinuities (nn take the 
form 

{pv^} = Q, [p + pvl] = (30) 

{B„} = 0, {vy} = 0, 

which corresponds to an ordinary hydrodynamic shock that propagates ac¬ 
cording to the conditions Byi = 0 and By 2 = 0 along the magnetic held. 
System (130|) is the transition solution between the oblique shocks in the limit 
01 —0 and the parallel shock 

As the angle 0i decreases, the discontinuity between the admissible mass 
huxes for the fast and the trans-Alfven shocks will also decrease. Conditions 
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ffTU]) in this case give 

B ^ 

{5^} = 0, {p'y^} = 0, [p + pv‘^] = -^, (31) 

} ^x2By2 1 P^x ■{ Vy } — BxBy2 ■ 

47r 

From the simultaneous solution of the last two equations, we have 

rri^ = P 2 BII 47r = (32) 

For this mass flux, Eq. mni has no unique solution. A nonzero Q\ can 
correspond to zero 02 - The tangential magnetic held component appears 
behind the shock front, corresponding to the switch-on shock (Son)- It is 
indicated in Fig. [2]a by the vertical segment at "^I^® switch-on 

shock can act as the transition one for the trans-Alfven and fast shocks in 

the limit 6*1 —>■ 0, but this, of course, requires that —)■ Otherwise, 

there will be the transition to the parallel shock according to fl30|l . 

To establish the form of the transition solution between the parallel shock 
and the contact discontinuity, we will set Ua, = 0 in (l30ll . After that we will 
obtain 

{B^} = {vy} = {p} = 0, (33) 

This system of equations corresponds to the contact discontinuity fl2T)) . or¬ 
thogonal to the magnetic held lines. It describes the transition discontinuity 
between the parallel shock in the limit Va, —)■ 0 and the contact discontinuity. 
Such a transition takes place at = 0 in Fig. Ela. 

When Bj; = 0 (Fig. inf>), becomes zero and all nonzero mass huxes 
locate in region IV in Fig. [2l To hnd the boundary conditions corresponding 
to them, let us substitute 5a; = 0 in (1T0|) . We obtain 

r 5 ^ 'i 

{pV^} = Q, + ^ -^1 = 0, (34) 

{ r>a:5y } = 0, {r>j^} = 0. 

These conditions characterize a compression shock propagating perpendicu¬ 
larly to the magnetic held. In the general case of a perpendicular shock {S±) 
we will hnd the boundary conditions by substituting 5a; = 0 in (fT|)-(]7|): 

{p'yx} = 0, + +-J 
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(36) 


{ VxBy ]■ 0 , Vy } 0 , 

{ V:,B^ } = 0, {i;^} = 0. 

Equations flM|l are then the boundary conditions for the transition discon¬ 
tinuity between the fast shock in the limit —)• 0 and the perpendicular 

shock with the magnetic field directed along the |/-axis. This transition can 
take place only at mass fluxes that satisfy inequality (11911 which takes the 
form > m\ at = 0 (Fig. |3]&). 

To determine the boundary conditions for the discontinuity at = 0 
(Fig. [3]c), let us substitute = 0 and = 0 in (II])-([7j). In this case, the 
magnetic held and the velocity held are parallel to the discontinuity surface 
and can undergo arbitrary jumps in magnitude and direction, while the jump 
in pressure is related to the jump in the magnetic held by the condition 


p + 


By^ + 
Svr 


= 0 , 


(36) 


This corresponds to the tangential discontinuity (T). The contact disconti¬ 
nuity, the slow shock, and the Alfven discontinuity can pass to it in the limit 
Sa, —>■ 0 under certain conditions. Let us hnd the corresponding transition 
solutions. Firstly, let us substitute 5a; = 0 in the boundary conditions for 
the contact discontinuity fl2Tl) . We will obtain the transition solution 

{By} = {vy} = {p} = Q. (37) 

It describes the tangential discontinuity fl36|l for the zero held component 
Bz in the absence of jumps {vy} and {By}. Secondly, the conditions for the 
oblique shocks (fTOj) at 5a; = 0 and = ^ are the transition solution 

{p+f}=0. (38) 


that corresponds to the plane tangential discontinuity fl36ll at 5^ = 0. Fi¬ 
nally, the boundary condition for the Alfven discontinuity fl29ll after the 
substitution of 5a; = 0 gives the transition solution 

{p} = 0, (39) 

that describes the tangential discontinuity fl36|l . with out any jump in density 

P- 
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Figure 4: Scheme of continuous transitions between MHD discontinuities. The dashed 
line encircl es the set of discon tinuities corresponding to the block of oblique shocks in the 
scheme bv ISvrovatskiil (195^. The shading inside the contour highl ights the slow (left) 
and fast (right) components of the scheme proposed bv ISomov (11994 . 


The angle 62 of the strongest trans-Alfven shock is dehned by Eq. (ET]). 
In view of (j2T|l . 62 —)■ —0i when 0. Thus, the trans-Alfven shocks 

degenerate into a special case of the Alfven discontinuity as the magnetic 
flux decreases. Of course, the density jump can also be set equal to zero 
for any type of flow. In this case, all parameters and will be 

equal to the flux (EH]), at which the Alfven discontinuity fl27ll takes place. In 
Fig. |3]c, it is denoted by m^oSAon- other mass fluxes, the differences in 
plasma characteristics on different sides of the discontinuity will disappear; 
the discontinuity will be absent as such. 

Let us combine the properties of discontinuous solutions systematized 
above into the scheme of permitted transitions shown in Fig. jH Here, the 
two-dimensional discontinuities are located in the middle row in order of in¬ 
creasing the mass flux and the three-dimensional discontinuities are located 
in the upper row. The one-dimensional parallel shock {S y) occupies the lower 
row. 

The hrst description of transition solutions f Svrovatskii . 1956[l contained 
only four types of discontinuous flows: a tangential discontinuity (T) and 
Alfven (A), oblique (S) and perpendicular {Sj_) shocks (Fig. |5|). The corre¬ 
sponding scheme of continuous transitions between discontinuous solutions 
of the equations of ideal MHD showed such transitions to be possible in prin¬ 
ciple, but it was dehnitely incomplete. Firstly, it did not have some of the 
discontinuous solutions, in particular, the parallel shock ( 511 ) and the con¬ 
tact discontinuity (C). Secondly, the block of oblique shocks (S) combined 


19 


























































Figure 5: 
19561 ^ 


Scheme of continuous transitions between MHD discontinuities l|Svrovatskii 



Figure 6: Scheme of continuous transitions between MHD discontinuities (jSomovl . Il994[ i 


several different discontinuities at once: fast (S'+) and slow (S'-) shocks, 
switch-on (Son) and switch-off (Sog), shocks, and trans-Alfven (Tr), shocks, 
the possibility of transitions between which requires a separate considera¬ 
tion. Subsequently, this picture of transitions was supplemented b ased on 
the c orrespondence between the shocks and small amplitude waves fjSomovi . 


1994f) . Although this approach allows the possible transitions (Fig. [6]) and 


even their conditions to be correctly specihed, it provides no description of 
the specihc form of transition solutions between the discontinuities under 
consideration. 

We group individual elements for the convenience of comparing the gen¬ 
eralized scheme of t r ansiti ons with those proposed previously. Syrovatskiis 
scheme fjSvrovatskiil . 119561) is consistent with Fig. H] if we combine the ele¬ 
ments Fi, Sos, Sl, Tr, S'on and 5'+ into one block of “oblique shocks” (S'), 
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while omitting the question of whether any transitions inside the block 
are possible and disregard the contact discon t inuity {C) and the parallel 
shock The scheme proposed in fISomovl . Il994l ). includes the parallel 

shock (S'II) and the separation of oblique shocks into the slow one (S'-), corre¬ 
sponding to condition flTSl) . and the fast one (S+) corresponding to condition 
flT^ . It is quite obvious that the scheme of transitions we propose is a proper 
and natural generalization of the two previous schemes. Our scheme contains 
not only evolutionary types of discontinuities but also non-evolutionary ones: 
the switch-on, Alfven, and trans-Alfven shocks. 


6. Jump in internal energy 

To determine the plasma heating efficiency, let us t urn to the boundary 
cond ition (jSj), which is the energy conservation law (ILedentsov fc Somovi . 
20131). Using Eq. ([2]), we will find the jump in internal energy from ([8]) 




{bV-(v-b)b4. 


(40) 


Using the mean velocities Vx, Vy and we will write the first term as 


= -Vx - Vy {vy} - {v^} 


We will express the jumps in tangential velocity components in terms of the 
jumps in tangential magnetic field components using Eqs. dl]), 02]), (E]) and 
dS]) as 

Now, the first term on the right side of Eq. fHOl) appears as 

,2 


1^1 _ (. .. VyBx r 75 -1 VzBx r 

“ IT J “ “ 4^ “ 4^ 

Let us substitute in the second term written as 

{vxP] = -—{vx} - —{p}, 
m mm 


(41) 


the pressure jump from Eq. (13), namely 

B 


B. 


{p} = -m{vx} - ^ {By} - ^ [B,]. 
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Here, as in Eq. fHD]) . we made use of condition (j2]). As the result, the second 
term in Eq. fl4Up takes on the form 


1 r T P r T ~ r T ^ 

— {Vxp} = ^ [By] + 

m m Airm 


VxBz f 


Let us now open the scalar product (v ■ B) in the third term: 

1 


dvrm 


{B'- 


y-B)B^} = 


(42) 


47rm 


{^xBy VyBx)By A (vxB^ ^ z B x} B z} ■ 


Then, let us apply conditions Eq. ([3]) and Eq. (|1]) to the derived equation: 


1 

dvrm 



(vB)H4 = 


"^xB^j VqjB. 


y-L-’x 


dvrm 


{B,}- 


^xBz VzBx 

Aura 


[B.]. 


(43) 


Thus, each of the three terms in the right side of Eq. fHU]) is expressed via 
individual jumps of the normal velocity components and tangential magnetic 
held components. Let us substitute Eqs. flTT]) - fl43|) in fl40|) 


{e I = — Kl + —— 

m Awm 


"^yBx 


{By}- 


VxBy 


dvrm 


VyBx r 1 

" -{By} + 


^xBz ^zBx f o T ^X.Bz 

+-^- {Bz} - 


^ zBx 


{B.}. 


dvrm dvrm 

This equation can be simplihed if we expand the means appearing in it: 





{Vx}{By} 

Ibvrm 


{B,} 


16nm 


{B,}. 


Factoring out —{vx}/'m = —{r}, we derive the hnal equation, which ex¬ 
presses the internal energy jump at the discontinuity through the jumps of 
the inverse density and tangential magnetic held components 


{e} 


-{r} 



{Byr + {BzY 
IQtt 


(44) 
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For two-dimensional discontinuities, it takes quite a simple form, 



( 45 ) 


Equation (m|) allows definitive conclusions regarding the change in plasma 
internal energy when crossing the discontinuity surface to be reached. First, 
the internal energy increases, because, according to Zemplens theorem, 
—{r} > 0 and p and {By}'^ are positive. Second, the change in internal 
energy consists of two parts: the thermodynamic and magnetic ones. The 
latter depends on the magnetic held conhguration and, hence, on the type 
of discontinuity. Let us express the tangential magnetic held components in 
Eq. fITS]) in terms of the corresponding inclination angles: 


{ e } = -{r} p - {r} ^ (tan 62 - tan ^i) ^ 


Then, we will take the thermodynamic part of the heating independent of 
the type of discontinuity as the zero point and will measure the jump in 
internal energy itself in units of —{r}B^/ IGtt. For this purpose, let us make 
the substitution 



We will obtain the equation 


{ e }^ = (tan 62 — tan 61 )'^ . 


(46) 


The magnetic part of the internal energy jump depends on the magnetic 
held conhguration and, as a result, on the discontinuity type. For the cal¬ 
culations shown in Fig. [21 the dependencies of the internal energy jump on 
the mass how passing through the discontinuity which were calculated using 
formula (H5|) . are shown in Fig. |6l Here, zero stands for the thermodynamic 
part of the jump, which does not depend on the mass how. It can be seen 
that, at the specihed plasma parameters, the maximum jump of the inter¬ 
nal energy is induced by the strongest trans-Alfven shock wave; moreover, 
its magnitude rapidly grows with an increase in the magnetic held incidence 
angle 61 . The correlations between the efficiency of plasma heating by other 
types of discontinuities depend on specihc conditions of the medium. For 
example, heating by slow shock waves can be both lower than heating with 
fast shock waves at smaller 61 , and higher than that at bigger di. In any 
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Figure 7: Jump in internal energy e versus mass flux for various angles 9i. 5° (solid lines), 
25°(dashed lines), and 45° (dotted lines). 

case, the heating depends on the shock strength. The larger the change in 
magnetic energy density, the higher the temperatures to which the plasma 
will be heated. 

The curve describing the jumps in internal energy at the strongest trans- 
Alfven and fast shocks is 

^ ^ ^ ’ 
' A 

It is indicated by the thin line in Fig. |6l 


7. The model of reconnection 


Markovskii fc SomovI fjl989l ) suggested a two-dimensional stationary 
model of magnetic recoi inection, t hat ge neral izes and combines the famous 
reconnection mod e ls by Petschelj f 19641) and Svrovatsl^ ( 1971 ). Later on. 


Bezrodnykh et al.l (120071) constructed a two-dimensional analytical model of 


reconnection in a plasma with a strong magnetic held that included a thin 
current layer {CL) and four discontinuous MHD hows of hnite length R at¬ 
tached to its endpoints. L is a halfwidth of the current layer. Figure [HI 
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displays the configuration of discontinuities inside which the electric currents 
flow perpendicular to the figure plane. The normal magnetic field component 
vanishes on the current layer and is equal to a given constant [3 on the MHD 
shock waves inclined to the x-axis at angle a. The field at infinity grows 
linearly with a proportionality coefficient 7 . The quantities L, R, a, /3, and 
7 are free parameters of this model used below. 

Magnetic field can be written as B = {Bx{x,y), By{x,y), 0), i.e., only 
two of its components, B^ and By, are nonzero and depend only on the x 
and y coordinates. It follows from the magnetic field potentiality that the 
function B = B^ — iBy is an analytic function of variable z = x + iy itl the 
exterior of the sections on the complex plane shown in Fig. [HI Thus, the 
model considered in this section is reduced to finding the analytic function 
B = Bx — iBy from the linear combination of B^ and By specified at its 
boundary. Since this problem is symmetric with respect to the x and y axes, 
it will suffice to consider it only in the first quadrant. Then the boundary 
condition for the function B is formulated as follows: By = Q on the current 
layer and surface (x = 0); Basina + By cos a = —(3 on the discontinuity 
surface; B^ = 0 on the su rface (y = 0, x > L) and B = i'jz + o(l) at 2 ; —)■ 00 . 


Bezrodnykh et ah! (l2007l) constructed the sought-for analytic function in an 


explicit form by the method of the conformal mapping of the domain outside 
the system of sections onto the half-plane. 

Figure [H] shows magnetic field lines calculated in the special case of re¬ 
connection {L = 1, R = 1, a = 45°, /5 = 1, 7 = 1) that can be taken as 
a basis for comparison with many other possible reconnection regimes. The 
figure demonstrates that the result of solving the problem in the main (cen¬ 
tral) part of the reconnection region is a current layer {CL). It is intersected 
by two symmetrically located magnetic fi eld lines; the points of intersection 
(for more details on their properties, see Svrovats^. 197ll; Somov . 2013bl ) 
separate the current layer areas with the field circulation relative to them 
having opposite signs. Thus, the Syrovatskii-type current layer that consists 
of a direct current and two attached reverse currents is actually located in 
the central part of the reconnection region. 

In specific astrophysical applications, in particular, in solar flares, the 
model of the so-called super-hot (Tg > 10 keV) turbulent-current layer 
f Somovl . 2013bh should be used to determine the physical parameters of this 
region. The advantage of an analytical model is the possibility to investigate 
more general patterns that do not depend on the detailed assumptions of a 
physical reconnection model. Consider some of the properties of discontinu- 
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Figure 8: The configuration of electric currents (thick straight segments) consists of a 
current layer (CL) and four segments (discontinui ty surfaces) of finite leng th R attached 
to its endpoints; L is the current layer half-width ( Bezrodnykh et al] . 120071) . 



Figure 9: Magnetic field lines (thin curves) with an indication of the values of the vector 
potential. The regions of direct and reverse currents inside the current layer (the thick 
horizontal straight segment) are seen. The magnetic field undergoes a.jump on the attached 


discontinuity surfaces (inclined straight segments) ( Bezrodnykh et al. . 20071) . 
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Figure 10: Distribution of the angles 0i and 82 at the attached shock. The variable I 
i s the distance measured from the current layer endpoint along the discontinuity surface 


(jBezrodnvkh et al.l . DOlll l. 


o us flows in the vicin i ty of a current layer predicted by the analytical model 
of Bezrodnykh et al. ( 2007 1. 

A magnetic held calculation in the approximation of a strong held near the 
current structure gives the angles of incidence and refraction of the held on 
the discontinuity surfaces from which the types of discontinuities are deter¬ 
mined. The characteristics of the plasma howing through the discontinuity 
change with the distance measured from the point of attachment to the cur¬ 
rent layer along the discontinuity surface. The set of these characteristics 
determines the type of discontinuity in each point of the attached surface. 
Moving from point to point along the surface of the discontinuity we will see 
the change of the flow conditions, and hence the change of the type of discon¬ 
tinuity. The different types of discontinuous MHD solutions will correspond 
to the different regions of the attached surface. The transition solutions will 
realize at the boundaries between t hese regions. _ 

For example, Fig. [TOl taken from iRezrodnvkh et al.l (1201111 shows how the 
magnetic held inclination angles change with coordinate at the discontinu¬ 
ity for one specihc calculation. The discontinuity surface can be arbitrarily 
divided into three regions in accordance with the types of hows. The trans- 
Alfven shock is immediately adjacent to the current layer. This is the domain 
of negative angles 61 in Fig. [101 As one recedes from the current layer, a fast 
shock is observed {62 > Oi) up to the point of intersection between the plots 
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of 6*1 and 6*2 at / = /q. The discontinuity ends with a slow shock ( 6*2 < 6 * 1 ) 
that gradually passes into a continuous flow ( 6*2 = 6 * 1 ) at / = 1. The absolute 
values of the angles 61 and 6*2 tend to the same value near the current layer 
(/ = 0 in Fig. fTOj) . Therefore, the Alfven discontinuity can be realized at this 
point. We also expect further continuous transitions at / = 0 in accordance 
with the scheme of transitions (Fig. Hj) but the accuracy of the Fig. [TUI does 
not allow us to say it for sure. 

As the plasma flow parameters change continuously, the type of shock 
must change through transition discontinuities. The hrst transition occurs 
between the trans-Alfven and fast shocks in the case we consider. As has 
been shown above, the switch-on wave fl3ip must act as a transition discon¬ 
tinuity. Indeed, as can be seen from Fig. [TD], the magnetic held at the point 
of the transition between the trans-Alfven and fast shocks is normal to the 
discontinuity surface ( 6*1 = 0 ) in the inhow and has a tangential component 
( 6*2 > 0) in the outhow. This transition cannot be made only by a gradual 
change in the mass hux through the discontinuity. The angle of incidence 
of the magnetic held should be simultaneously reduced to reduce the dis¬ 
continuity between the admissible mass huxes for the fast and trans-Alfven 
shocks. In addition, both trans-Alfven and fast shocks satisfy inequality (IT^ . 
Hence, for the parameters of the medium to change smoothly, the mass hux 
through the switch-on wave must also satisfy this inequality: 

2 ^ 2 I 2 

^oS > + m_L , 


or 


B: 


> 


B^ + B- 


47rr2 dvrr 

Simplifying this relation, we obtain a constraint on the possible inclination 
of the magnetic held behind the discontinuity plane: 


tg2 02<^i^. (48) 

Pi 

The second transition between the fast and slow shocks at Z = Zq cannot 
be made by a continuous change in the mass hux through the discontinuity. 
The intersection of the curves in Fig. [10] suggests that there is no jump in 
magnetic held strength at the point separating the fast and slow shocks. Only 
the contact discontinuity can correspond to such a held structure. However, 
as has been shown above, the ratio of the tangents of the angles tan 6 * 2 /tan 0 i 





in the fast shock asymptotically tends only to P 2 /Pi from above and can 
take on unity typical of a contact discontinuity only at {p} = 0. This, along 
with the boundary conditions fl2T|) . gives the absence of a discontinuity at 
the point of contact between the fast and slow shocks. The jumps in density, 
velocity, magnetic held, and pressure at this point must be zero. 

Thus, the discontinuity surface turns out to be physically separated into 
two regions: the inner part consists of trans-Alfven and fast shocks, while 
the outer part is a slow shock. When changing the initial model parameters, 
we can obtain a reconnection regime in which the discontinuity ends with 
a fast shock, while the outer part of the discontinuity is absent altogether. 
This suggests that the inner part of the discontinuity is attributable to the 
reconnection process itself and is closely related to the presence of a reverse 
current at the endpoints of the current layer, which was clearly shown by 
Bezrodnykh et ah fl2007f) . while the outer part depends strongly on the fac¬ 


tors affecting the overall topology of current layers: the presence or absence 
of a “magnetic obstacle”, non-uniformity of the plasma distribution in the 
reconnection region. 

We can trace some analogies of our conclusions wit h the results of present- 


day numerical MHD simulations of fast reconnection (IShimizu fc Ugail. 12003 


Ueai et ah . 2005 : Ugai, 20081: Zenitani fc T. Mivoshi . 2011 1. As a result of 


reconnection, the plasma being ejected from a current layer is gathered into 
the so-called “plasmoid” separated from the surrounding plasma by a system 
of slow shocks. The structure and intensity of the latter depend on the plas¬ 
moid sizes and density. In addition, the chara c terist ics of the outer part of 
the discontinuity surface in iBezrodnvkh et al.l (1201 ll) depend on the model 
geometry. Closer to the current layer, in the region of reverse currents, the 
numerical models give a complex system of shocks. Since the derived trans- 
Alfven waves are non-evo lutionary, the structure of the inner part of the 
discontinuity described by IBezrodnvkh et al.l (120071 1201 ill must also become 
more complicated and can have similarities to the results of numerical simu¬ 
lations. On the whole, the interrelationship between the process of magnetic 
reconnection and the formation of a system of accompanying discontinuous 
flows requires a further comprehensive study. 

Equation (H^ allows us to draw some dehnite conclusions concerning vari¬ 
ations in the internal plasma energy upon the passage though the disconti¬ 
nuity surface. First, the internal energy increases since, according to the 
Zemplen’s theorem, —{r} > 0, and p and {ByY are positive quantities. 
Second, the internal energy change is made of two parts, namely, a thermo- 
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Figure 11: RHESSl images in different energy bands at 23:11:00 UT. The three contours 
(80% of the peak flux in each image) on the solar disk indicate the top of loops in the energy 
bands (contour line shade from light to dark): 6-8, 10-12, and 16-20 keV. The contours 
(80% of the peak flux of the coronal source in each image) above the limb are for the (from 
light to dark) 10-12, 12-14, and 14-16 keV bands. The asterisk marks the centroid of each 
source. The crosses mark the two footpoints of the X-ray loop ( Sui fc Holman . 2003ll . 


dynamic part, which is defined by the plasma pressnre and a magnetic one, 
which is related to variations in the magnetic field strncture in the vicinity 
of the discontinnity surface. 

The maximum heating at discontinuities can be expected in a magnetic 
field with a sharply varying geometry. These conditions are met in the re¬ 
gion of magnetic reconnection. Similar discontinuity structures can be intro¬ 
duced as independent elements into analytical reconnection model. Earlier 
we showed that the waves that are near the face ends of the current layer 
(where reverse current are generated) are trans-Alfven shock waves. Moving 
away from the current layer, the jump of the magnetic field intensity at the 
discontinuity decreases and the density jump vanishes. Thus, the best condi¬ 
tions for plasma heating occur near the region of reverse currents. However, 
this conclusion is purely qualit ative. Since the input parameters of the model 
Bezrodnykh et ahl fl2007l. I 2 OIII) are some abstract quantities, we can not pro¬ 
duce any numerical estimates about the real plasma heating. This process 
requires further study. Nonetheless, heating by shock waves can help form 
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a s uper-hot plasma obse rved with modern X-ray space observatories. (Fig. 
[TTl Sui fc Holman . 2003 '). 


8. CONCLUSIONS 

The correspondence between the standard classihcation of discontinuous 
flows in ideal MHD and the characteristic parameter of plasma flow - the 
value of the mass flow through the discontinuity - is set up. On this basis, 
all of allowed transition solutions are found. As a result, the generalized 
scheme of the continuous transitions between MHD discontinuities is built. It 
contains the discontinuities that were not represented in the earlier schemes, 
such as the contact discontinuity, swich-on and swich-off shocks. Some types 
of discontinuous flows, for example the trans-Alven shock waves, are non¬ 
evolutionary. They are also included in the generalized scheme of transitions. 

Within a framework of the simplihed analytical model of the magnetic 
reconnection various areas of connected to the current layer discontinuity 
surface are identihed with different types of MHD shock waves. In particu¬ 
lar, regions of the trans-Alven shock waves are located near the ends of the 
current layer (in the presence of reverse currents). Division of the attached to 
the current layer discontinuity surfaces into two areas is found. That division 
is result of the different origin of discontinuities. The quasi-st at ionary inner 
region is due to the reverse currents in the current layer, and the outer region 
is caused mainly by the boundary conditions of the magnetic reconnection 
process and by the speed of reconnection. 

We have considered dependence of the plasma heating as the thermody¬ 
namic parameters of plasma as the type of MHD discontinuity. The larger the 
jumps in plasma density and magnetic energy density at the discontinuity, 
the stronger the heating. Such conditions occur near a region of the reverse 
currents in the reconnection process. We believe that this result will be use¬ 
ful for quantitative explaining the temperature distributions of the super-hot 
plasma in solar flares observed by modern space X-ray observatories. 
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